Switching and defect dynamics in multistable liquid crystal devices 
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We investigate the switching dynamics of multistable nematic liquid crystal devices. In particular we identify a 
remarkably simple 2-dimensional (2D) device which exploits hybrid alignment at the surfaces to yield a bistable 
response. We also consider a 3-dimensional (3D) tristable nematic device with patterned anchoring, recently 
implemented in practice, and discuss how the director and disclination patterns change during switching. 
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Liquid crystals (LCs) are interesting soft materials 
with widespread technological applications, for instance 
liquid crystal displays (LCDs). The archetype LCD de- 
vice is the "twisted nematic display", commonly em- 
ployed in flat panel monitors^. This device is easy to 
build in practice, but its functioning is not optimal. 
Firstly, it requires the constant application of an electric 
field for the whole time the device is to be in the "on" 
state. Secondly, standard twisted nematic LCD devices 
suffer from the so-called "viewing angle problem" 2 , which 
refers to the limited range of angles from which the dis- 
play is readable. Both these shortcomings can be avoided 
at least conceptually. On one hand, the viewing angle 
range can be dramatically enhanced by building a "multi- 
domain" device, with the LC anchoring patterned on the 
device surface: this leads to the simultaneous presence 
of regions of left and right-handed twist, whose viewing 
angle characteristics differ and complement each other. 
Modelling plays an important role in optimising the deli- 
cate balance between the gain in viewing angle range and 
the loss in resolution due to the creation of defects be- 
tween regions with different director orientation^-. On 
the other hand, bistable or multistable devices can re- 
tain memory of two or more distinct states, even after 
the electric field is switched off. These solved the prob- 
lem of having to constantly use a large field to keep the 
device in the "on" state. Multistable devices are also in- 
teresting theoretically as they require the presence of a 
series of met ast able states, with comparable free energy 
yet different defect topology and optical properties. How 
to realise this in practice is usually non-trivial. Two well 
known multistable devices proposed in the physical liter- 
ature are the zenithal-bistable device^, and the tri-stable 
device experimentally demonstrated in 8 . In both cases, 
defect dynamics is likely to be important^. 

Here we describe the physics behind the switching dy- 
namics of two multistable and multidomain nematic de- 
vices. Our main result is that it is possible to design 
a remarkably simple 2-domain hybridly aligned nematic 
(HAN) cell (Fig. la). For a wide range of parameters we 
find two metastable states with competing free energy: 
one in which the LC pattern is made up by stripes with 
alternating splay-bend direction with defects in between, 
and another one in which the device becomes effectively 
a single domain with either clockwise or anticlockwise 
deformation. The dynamical schedule with which the 
electric field is switched on and off is crucial to ensure 



bistability. We also investigate a 3D tristable device sim- 
ilar to the one built in 8 - (Fig. lb), and our results suggest 
that the disclination dynamics involved in the switching 
is strikingly non-trivial. Finally, we quantify the effect 
of backflow in the switching of the first of these devices 
and we show that it can change the range of voltages for 
which multist ability occurs. 

The physics of the LC device may be described by a 
Landau-de Gennes free energy dependent on a tensor or- 
der parameter, Q a p, whose largest eigenvalue (denoted 
by 2/3q) gives the magnitude of the local order. The free 
energy density / is a sum of three terms. The first is a 
bulk contribution, 
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where Ao is a constant and 7 controls the magnitude 
of the ordering (Greek indices denote Cartesian compo- 
nents and summation over repeated indices is implied). 
A second term K/2 (d 1 Q a f3) 2 describes the free energy 
cost of distortions^ where K is an elastic constant. The 
last term, —^£-E a Q a pEp, describes the interaction with 
the external electric field E where e a > is the dielec- 
tric anisotropy of the material. The equation of motion 
for Q is D t Q = TH, where T is a collective rotational 
diffusion constant and D t is the material derivative for 
rod- like molecules 9 . The molecular field H ensures that 
Q evolves towards a minimum of the free energy^. The 
velocity field obeys a Navier- Stokes equation with a stress 
tensor generalised to describe LC hydro dynamics^i^. To 
solve the equations of motion, we used a hybrid lattice 
Boltzmann algorithm, see e.giii. In what follows we ex- 
press parameters in simulation units. To convert these 
to real units, one may assume to model a ~ /im thick 
device and a flow- aligning 1,9 LC with rotational viscos- 
ity equal to 1 poise - a typical value for e.g. 5CB. In 
this way, one force, time, and space unit may be mapped 
to 62.5 pN, 0.7 /is, and 0.025 /im respectively. The elec- 
tric field may be quantified via the dimensionless number 
e 2 = 27e a E^/(327rAoj). The simulation domains along 
the three coordinate axes are named L x , L y and L z re- 
spectively (there are periodic boundaries along x and y). 

We begin by considering the two-domain device shown 
in Fig. la, with mixed anchoring, homeotropic - along z 
- at the top (z = L) and homogeneous - along y - at the 
bottom (z = 0) of the cell. A small pretilt is patterned to 
create a defect in the middle of the cell, where a region 
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FIG. 1. Geometry of the multidomain devices in this work, 
(a) 2-domain cell, with homeotropic anchoring at the top plate 
and homogeneous anchoring with pretilt patterning at the 
bottom plate, (b) Surface patterning used for the top and 
bottom boundary conditions of the 3-domain tristable device. 



of clockwise splay-bend director rotation meets a mir- 
ror image with anticlockwise deformation. We call this 
zero-field free energy minimum the "two-domain" state. 
The geometry of this device is 2D, as there is transla- 
tional invar iance along x. This device is a simple variant 
of the commonly used HAN cell, and it may be built 
by either chemical or topographical patterning. Our de- 
vice contains coexisting homeotropic and homogeneous 
domains, which might lead to a switching in either di- 
rection, through the growth of one or the other domain: 
hence bistability is possible. Fig. 2a-c shows the evolu- 
tion of the director field in response to an electric field 
which is switched on and off first along y, then along z. 
After application of the field along the disclinations 
separating neighbouring clockwise and anticlockwise do- 
mains annihilate and transform the two-domain device 
into a single domain one. This defect-free state is sta- 
ble after the removal of the field, as the pretilt is not 
enough to bend the director field in the bulk. On the 
other hand, when a sufficiently strong field is applied 
along z, defects reappear close to the surface, between 
regions with different pretilt. After switching off, elas- 
ticity alone is not enough to melt these defects away - 
rather they migrate back to the centre of the sample and 
reconstitute the initial two-domain structure. The time 
evolution of the free energy during this electric field cy- 
cle is shown in Fig. 2 (top panel, top set of curves) - this 
shows that both the single domain and the two-domain 
states are metastable and close in free energy, the main 
requirement for a bistable device. 

Changing the device size and the elastic constant af- 
fects both the switching on and off times and the bista- 
bility. E.g. relaxation times after switching off of the 
horizontal electric field, (see caption Fig. 2), scale like 
L 2 Z /KT (data not shown), as for common nematic de- 
vices^. The device is bistable only for small size (e.g. 
L z < 60 for K = 0.08 and other parameters as in the 
caption of Fig. 2) - for larger devices the application of 
a field perpendicular to the boundary fails to turn the 
single domain state back into the two-domain one. Fur- 
themore, bistability occurs for intermediate values of K 
- between about 0.06 and 0.09 for L z = 50, and other 
parameters as in Fig. 2 (Fig. 2, curves for K = 0.04 with 




FIG. 2. Free energy evolution (top panel) of the 2D device. 
Parameters were Aq = 1, T = 0.33775, 7 = 3, L y = L z = 40, 
e 2 - 0.06, K = 0.08 (open circles) and K = 0.04 (filled 
circles). All quantities are given in simulation units. Filled 
squares correspond to a simulation with backflow. The direc- 
tor field configuration of the corresponding stable states are 
reported in the bottom panel. The device is only bistable for 
K = 0.08. The step function in the top panel is non-zero 
when the field is applied. From t = up to t = 30000 the 
system relaxes from the configuration depicted in Fig. 1 to 
a stable state, ((a) in the bottom panel). An electric field 
along y is switched on at t = 30000 and off at t = 40000; 
the system relaxes to a stable state, ((b) in the bottom panel, 
snapshot at t = 70000) characterized by a new free energy 
minimum. At t — 70000 an electric field along z is switched 
on, and off at t — 80000. The relaxation process brings the 
system back to state (a) (with the defect in a slightly differ- 
ent position) . The relaxation curves for the two cases with no 
backflow (open and filled circles) were fitted via a stretched 

exponential, A + Be~^^ T ^ , (solid thick lines, see inset for 
relaxation to state (b) with K = 0.08, where the dashed line 
is instead a fit to a simple exponential). For K — 0.08, r and 
P are - 3200 and 0.72 (a), 540 and 0.52 (b); 80 and 0.39 (c). 
For K = 0.04, r and f3 are - 5140 and 0.51 (a), 3460 and 1.50 
(b); 360 and 0.56 (c). 



only one stable state). This is because when K is too 
small, there is not enough elastic driving force to lead to 
defect annihilation during the first half of the voltage cy- 
cle, whereas when K is too large the single domain state 
becomes too low in free energy with respect to the two- 
domain one, and the dynamics always gets stuck there. 

Backflow quantitatively affects the extent of the 
bistable regime. When it is neglected, the electric field 
range for which bistability is observed is overestimated. 
E.g., with no backflow, K = 0.08 and other parameters 
as in Fig. 2, bistability requires e 2 > 0.03, whereas with 
backflow e 2 > 0.05 (filled squares in Fig. 2, top). This 
is because hydrodynamics decreases the chance of get- 
ting stuck into the two competing metastable states with 
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FIG. 3. Evolution of the director field in the midplane (left 
column) and disclinations (right column) in the 3D tristable 
device. Snapshots corresponds to time equal to (a) 100, (b) 
3700, (c) 4600 and (d) 6600 simulation units. Parameters are: 
L x 90, L y 52, L z = 40, K = 0.08, e 2 0.0018, and other 
parameters as in Fig. 2. The field is applied at t = and at 
an angle of 60° with respect to x. 



different free energies. 

The idea of using domains of different director orienta- 
tions as seeds for multistable field-induced steady states 
has been used for the tristable device in£. The relevant 
geometry is shown in Fig. lb - it consists of a hexago- 
nal tiling of one or both boundary planes. In order to 
tile the whole boundary, three director orientations are 
needed. Via an electric field in the xy plane, the authors 
o& managed to stabilise each of the orientations used for 
the surface patterning. But what exactly is the switch- 
ing dynamics of this multistable device? What is the LC 
texture corresponding to the zero-field case? And what is 
the role of defect motion for the device switching? Such 
questions cannot be easily answered experimentally - yet 
are important to better control the device performance. 
Our numerical approach is ideal to fill in this gap and 
follow the director dynamics during switching (Fig. 3). 
For e 2 = 0, the surface patterning forces the formation of 
columnar straight disclinations - of both integer and half- 
integer topological charge (Fig. 3a). This 3-domain state 
is metastable. When an electric field is applied, along 
one of the anchoring directions, the line of asters, which 



marks the point at which six different domains meet, un- 
zips into a loop of half-integer topological charge. This 
dynamics is linked to a complex reshaping of the direc- 
tor field in the bulk (see Figs. 3b and 3c). Later on, the 
disclinations loop merge and reconstruct into an array 
of disclinations, which are now parallel to the boundary 
planes and pinned to the surface (Fig. 3d). This electric- 
field induced state is stabilised by the bulk elastic term 
even after switching off. Two important control param- 
eters determining the physics of our tristable device are 
the ratio between domain and sample size, d/L, and the 
"interface number" dj ^/ Aq/K proportional to the ratio 
between the domain and defect core size. Briefly, addi- 
tional simulations show that, if the sample is too thin, 
boundaries always restore the same state after switch- 
ing off the inplane field, and multistability is lost. On 
the other hand, increasing the "interface number" would 
stabilise a disclination network in the bulk of the device 
- however this would only become relevant for d in the 
sub-/im range, difficult to reach in practice. 

Thus, we studied the switching dynamics of multi- 
stable nematic LCs. These are important for applica- 
tions, as they enable a device to function without the 
constant application of an electric field. We have shown 
that a 2D HAN cell, with homeotropic anchoring at one 
boundary and planar patterned surface anchoring at the 
other boundary provides a strikingly simple example of 
a bistable device. The simplicity of the geometry con- 
sidered and of the electric field cycle should render this 
device easy to be built. We also elucidated the role of 
defect dynamics in the switching of a tristable nematic 
device akin to that of&, which was again stabilised by sur- 
face patterning. The spectacular defect restructuring in 
Fig. 3 should leave an observable signature in the optical 
patterns which can be recorded experimentally. 
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